The 155-day periodicity of the sunspot area fluctuations 
in the solar cycle 16 is an alias 
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Abstract The short-term periodicities of the daily 
sunspot area fluctuations from August 1923 to Octo- 
ber 1933 are discussed. For these data the correlative 
analysis indicates negative correlation for the periodic- 
ity of about 155 days, but the power spectrum analysis 
indicates a statistically significant peak in this time in- 
terval. A new method of the diagnosis of an echo-effect 
in spectrum is proposed and it is stated that the 155- 
day periodicity is a harmonic of the periodicities from 
the interval of [400, 500] days. 

The autocorrelation functions for the daily sunspot 
area fluctuations and for the fluctuations of the one 
rotation time interval in the northern hemisphere, sep- 
arately for the whole solar cycle 16 and for the max- 
imum activity period of this cycle do not show differ- 
ences, especially in the interval of [57, 173] days. It 
proves against the thesis of the existence of strong pos- 
itive fluctuations of the about 155-day interval in the 
maximum activity period of the solar cycle 16 in the 
northern hemisphere. However, a similar analysis for 
data from the southern hemisphere indicates that there 
is the periodicity of about 155 days in sunspot area data 
in the maximum activity period of the cycle 16 only. 

Keywords Sun: sunspot area fluctuations - Sun: mid- 
term periodicities 



1 Introduction 

For about 20 years the problem of properties of short- 
term changes of solar activity has been considered 
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extensively. Many investigators studied the short- 
term periodicities of the various indices of solar activ- 
ity. Several periodicities were detected, but the pe- 
riodicities about 155 days and from the interval of 
[470,620] days ([1.3,1.7] years) are me ntioned most 
often. First of them was discovered by iRieger et al 



( 19841) in the occurence rate of gamma-ray flares de 
tected by the gamma-ray spectrometer aboard the So- 
lar Maximum Mission (SMM). This periodicity was 
confirmed fo r other solar flares data and for the same 
time period feogart fe Bailll985l : feai fe Sturrocklll987 : 



Kile fe Oliver! 1 199 II ). It was al so found in proton flares 



during solar cycles 19 and 20 (|Bai fc Oliver! 1 1 9901 ) . but 
it was no t found in the solar flares data during solar 
cycle s 22 (jKile fc Cliverlll99ll: lBa3ll992t lOzgiic fc Atac 
19891) . 



Several autors co nfirme d abov e results for the daily 
sunspot area data. Lean ( 19901 ) studied the sunspot 
data from 1874-1984. She found the 155-day periodic- 
ity in data records from 31 years. This periodicity is al- 
ways characteristic for one of the solar hemispheres (the 
southern hemisphere for cycles 12-15 and the north- 
ern hemisphere for cycles 16-21). Moreover, it is only 
present during epochs of maximum activity (in episodes 
of 1-3 years). 

Similar investigations w ere carried out by 



Carbonell fc Ballesterl (jl992l ). They applied the same 



power spectrum method as Lean, but the daily sunspot 
area data (cycles 12-21) were divided into 10 shorter 
time scries. The periodicities were searched for the fre- 
quency interval 57-115 nHz (100-200 days) and for each 
of 10 time series. The authors showed that the period- 
icity between 150-160 days is statistically significant 
during all cycles from 16 to 21. The considered peaks 
were remained unaltered after removing the 11-year cy- 
cle and applying the power sp e ctrum analysis. 



Oliver. Ballester fc Baudinl (|1998l ) used the wavelet 



technique for the daily sunspot areas between 1874 and 
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1993. They determined the epochs of appearance of this 
periodicity and concluded that it presents around the 
maximum activity period in cycles 16 to 21. Moreover, 
the power of this periodicity started growing at cycle 
19, decreased in cycles 20 and 21 and disappered after 
cycle 21. 

Similar analyses were pr esented by 
Prabhakaran Navar et al.l (|2002 ) , but for sunspot num- 
ber, solar wind plasma, interplanetary magnetic field 
and geomagnetic activity index A p . During 1964- 
2000 the sunspot number wavelet power of periods 
less than one year shows a cyclic evolution with 
the phase of the solar cycle. The 154-day period is 
prominent and its strenth is stronger around the 
1982-1984 interval in almost all solar wind param- 
eters. The existence of the 156- day periodicity in 
sunsp ot data were confirmed by iKrivova fe Solanki 
(|2002l ). They considered the possible relation be- 



tween the 475-day (1.3-year) and 156-day periodici- 
ties. The 475-day (1.3-year) periodicity was also de- 
tected in variations of the interplanetary magnetic 
field, geomagnetic ac tivity helioseismic data and in 
the s o lar wind speed (IPaularena. Szabo fc Richardson 



1995t ISzabo. Lepping fc Kind Il995l: iLockwoodl 12001 



Rich ardson et al.l I1994I ). iPrabhakaran Navar et a!. 
( 20021 ) concluded that the region of larger wavelet 
power shifts from 475-day (1.3-year) period to 620- 
day (1.7-year) period and then back to 475-day (1.3- 
year). The periodicities from the interval [475,620] 
days ([1.3,1.71 year s) have been considered from 1968. 
Yacob fc Bhargaval (|1968l ) mentioned a 16.3-month 
(490-day) periodicity in the sunspot numbe rs and in 
the geomagnetic data. Bai fc Sturrockl (1987) analysed 
the occurrence rate of major flares during solar cycles 
19. They found a 18-month (540-day ) periodicity in 



flare r ate of the norhern hemisphere. Ilchimoto et al 



dl985l) confirmed this result for the H a flare data for 
solar cycles 20 and 21 and found a peak in the power 



spectra near 510-540 days. lAkioka et al.1 (| 19871 ) found 



a 17-month (510-day) periodicity of sunspot groups and 
their areas from 1969 to 1986. These authors concluded 
that the length of this period is variable and the reason 
of this periodicity is still not un derstood. 

Carbonell fc Ballesten Jl992l) an d 
Oliver. Ballester fc Baudin ( 19981 ) obtained statisti- 



cally significant peaks of power at around 158 days for 
daily sunspot data from 1923-1933 (cycle 16). In this 
paper the problem of the existence of this periodicity 
for sunspot data from cycle 16 is considered. The daily 
sunspot areas, the mean sunspot areas per Carrington 
rotation, the monthly sunspot numbers and their fluc- 
tuations, which are obtained after removing the 11-year 
cycle are analysed. In Section 2 the properties of the 



power spectrum methods are described. In Section 3 
a new approach to the problem of aliases in the power 
spectrum analysis is presented. In Section 4 numeri- 
cal results of the new method of the diagnosis of an 
echo-effect for sunspot area data are discussed. In Sec- 
tion 5 the problem of the existence of the periodicity 
of about 155 days during the maximum activity period 
for sunspot data from the whole solar disk and from 
each solar hemisphere separately is considered. 



2 Methods of periodicity analysis 

To find periodicities in a given time series the power 
spectrum analysis is applied. In this paper two meth- 
ods are used: The Fast Fourier Transformation al- 
gorithm with the Hamming window function (FFT) 
and the Blackman-Tukey ( BT) power spectrum method 
(jBlackman fc Tukevllll)58t ). 

The BT method is used for the diagnosis of the rea- 
sons of the existence of peaks, which are obtained by the 
FFT method. The BT method consists in the smooth- 
ing of a cosine transform of an autocorrelation function 
using a 3-point weighting average. Such an estimator is 
consistent and unbiased. Moreover, the peaks are mi- 
correlated and their sum is a variance of a considered 
time series. The main disadvantage of this method is 
a weak resolution of the periodogram points, partic- 
ularly for low frequences. For example, if the auto- 
correlation function is evaluated for i = 1,...,1000, 
then the distribution points in the time domain are: 
1000, 500, 400, 333, 250, 167, 154, . . . Thus, it is obvious 
that this method should not be used for detecting low 
frequency periodicities with a fairly good resolution. 
However, because of an application of the autocorre- 
lation function, the BT method can be used to ver- 
ify a 'reality' of peaks which are computed using a 
method giving the better resolution (for example the 
FFT method). 

It is valuable to remember that the power spectrum 
methods should be applied very carefully. The difficul- 
ties in the interpretation of significant peaks could be 
caused by at least four effects: a sampling of a contin- 
ues function, an echo-effect, a contribution of long-term 
periodicities and a random noise. 

First effect exists because periodicities, which are 
shorter than the sampling interval, may mix with longer 
periodicities. In result, this effect can be reduced by an 
decrease of the sampling interval between observations. 

The echo-effect occurs when there is a latent har- 
monic of frequency a in the time series, giving a spec- 
tral peak at a, and also periodic terms of frequency 
ia, i«, etc. This may be detected by the autocorrela- 
tion function for time series with a large variance. 
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Time series often contain long-term periodicities, 
that influence short-term peaks. They could rise pe- 
riodogram's peaks at lower frequencies. However, it is 
also easy to notice the influence of the long-term peri- 
odicities on short-term peaks in the graphs of the au- 
tocorrelation functions. This effect is observed for the 
time series of solar activity indexes which are limited 
by the 11-year cycle. 

To find statistically significant periodicities it is rea- 
sonable to use the autocorrelation function and the 
power spectrum method with a high resolution. In the 
case of a stationary time series they give similar re- 
sults. Moreover, for a stationary time series with the 
mean zero the Fourier transform is equivalent to the co 



sine transform of an autocorrelation function (Anderson 



19711 ). Thus, after a comparison of a periodogram with 
an appropriate autocorrelation function one can detect 
peaks which are in the graph of the first function and do 
not exist in the graph of the second function. The rea- 
sons of their existence could be explained by the long- 
term periodicities and the echo-effect. Below method 
enables one to detect these effects. 



3 A method of the diagnosis of an echo-effect 
in the power spectrum 

The method of the diagnosis of an echo-effect in the 
power spectrum (DE) consists in an analysis of a pe- 
riodogram of a given time series computed using the 
BT method. The BT method bases on the cosine 
transform of the autocorrelation function which creates 
peaks which are in the periodogram, but not in the 
autocorrelation function. 

The DE method is used for peaks which are com- 
puted by the FFT method (with high resolution) and 
are statistically significant. The time series of sunspot 
activity indexes with the spacing interval one rotation 
or one month contain a Markov-type persistence, which 
means a tendency for the successive values of the time 
scries to 'remember' their antecendent values. Thus, 
I use a confidence level bas ing on the 'red noise' of 
Markov ( Mitchell et al. 19661 ) for the choice of the sig- 
nificant peaks of the periodogram computed by the 
FFT method. When a time series does not contain 
the Markov-type persistence I apply the Fisher test 
and the Kolmog orov-Smirnov test at the significance 
level a = 0.05 ( Brockwell &: Davis] 1991 ) to verify a 
statistically significance of periodograms peaks. The 



Fisher test checks the null hypothesis that the time 
series is white noise agains the alternative hypothesis 
that the time series contains an added deterministic 
periodic component of unspecified frequency. Because 
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Fig. 1 Two upper curves show the periodogram of the 
time series {S n (ti) — S„(U)} (solid line) and the 95% con- 
fidence level basing on the'red noise' (dotted line). The pe- 
riodogram values are presented on the left axis. The lower 
curve illustrates the autocorrelation function of the same 
time series (solid line) . The dotted lines represent two stan- 
dard errors of the autocorrelation function. The dashed 
horizontal line shows the zero level. The autocorrelation 
values are shown in the right axis. 
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Fig. 2 Same as in Fig. 1, but for the time series {S 3 (ti) — 
S s (ti)}. Because the statistical tests indicate that the time 
series is a white noise the confidence level is not marked. 
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the Fisher test tends to be severe in rejecting peaks as 
insignificant the Kolmogorov-Smirnov test is also used. 

The DE method analyses 'raw' estimators of the 
power spectrum. They are given as follows 



d = {t : c T > A cos( — ) > A (21 — 1)— < 
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r<(2Z + l)^}, for l = l,...,k-l, (5) 
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(1) C fc = {r : c T > A cos(— ) > A (2k - 1)— < 



for k=l,...,M, 
where 
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A 



for t — 1 , . . . , M, 

N is the length of the time series {Xi} and X is the 
mean value. 

The first term of the estimator S k is constant. The 
second term takes two values (depending on odd or even 
k) which are not significant because cm ~ for large M. 
Thus, the third term of (1) should be analysed. Looking 
for intervals of r for which ^2 T c T cos(^f) has the same 
sign and different signs one can find such parts of the 
function c T which create the value Sk- 

Let the set of values of the independent variable of 
the autocorrelation function be called / = {r : r = 
1, . . . , M} and it can be divided into the sums of disjoint 
sets: 



i = \jA t u \jBj u \jc h 

where 

A = {t : c T cos(^f ) < A < r < f£ }, 



.7rfer, (2i-l)M 
Ai = {t:c t cos( — ) < A i — < r < 



(2) 



(2i + l)M 
2k 



} , for i = 1, . . . ,k — 1, 



A k = {t: c T cos(^) < A (2k — < r < M} 



B = {t: c t < A cos(^) < A < r < 



M 



Bj = {t : c T < A cos(— — ) < A (2j - 1)— < 



AI 



2k 



M. 



r<(2i + l)— }, for i = l,...,fc-l, (3) 
S fc = {r : c T < OAcos(^) < 0A(2/c-l)|f < r < M}, 

7T k 1" ^\'jT 

Co = {t : c T > Acos( — ) > 0A0 < r < — }, (4) 
M Ik 



< AI}. 
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Well, the set Ci contains all integer values of r from 
the interval of ((2/ - l)|f , (2/ + ) for which the au- 
tocorrelation function and the cosinus function with the 
period 2j£ arc positive. The index / indicates succes- 
sive parts of the cosinus function for which the cosinuscs 
of successive values of r have the same sign. However, 
sometimes the set Ci can be empty. For example, for 
k = 13 and M = 1000 the set C\ should contain all 
r e [38, 115] for which c T > and cos(^f) > 0, but 
for such values of r the values of cos( :!! jf) are negative. 
Thus, the set C\ is empty. 

Let us take into consideration all sets {Ai}, {Bj} 
and {C{} which are not empty. Because numbcrings 
and power of these sets depend on the form of the au- 
tocorrelation function of the given time series, it is im- 
possible to establish them arbitrary. Thus, the sets of 
appropriate indexes of the sets {A t }, {Bj} and {C;} are 
called I a, Ib and Ic respectively. For example the set 
Ic contains all I from the set {1, . . . , k} for which the 
sets C\ are not empty. 

To separate quantitatively in the estimator Sk the 
positive contributions which are originated by the cases 
described by the formula (5) from the cases which are 
described by the formula (3) the following indexes are 
introduced: 
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9 k (l) = c T cos(-^-) for each I € /< 
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Fig. 3 Same as in Fig. 1, but for the time series {R(ti) — 



ws k (i) = X! Ct cos ^ _ /w - ^ for each :) ' e Ib ' 

reBj 

WS k = ^ W s kti) + WS k( 1 )' 

jei B i&ic 

taking for the empty sets {Bj} and {Ci} the indices 
ds^(l) and ds^(l) equal zero. 

The index ds^(l) describes a percentage of the con- 
tribution of the case when c T and cos(^f) are positive 
to the positive part of the third term of the sum (1). 
The index ds^ (I) describes a similar contribution, but 
for the case when the both c T and cos( z ^ : ) are simul- 
taneously negative. Thanks to these one can decide 
which the positive or the negative values of the auto- 
correlation function have a larger contribution to the 
positive values of the estimator Sk- When the differ- 
ence ws^(l) — ws^(j) is positive, the statement 'the 
fc-th peak really exists' can not be rejected. Thus, the 
following formula should be satisfied: 

2ds+ - 100 > 0. (7) 

Because the k-th peak could exist as a result of the 
echo-effect, it is necessary to verify the second condi- 
tion: 
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Fig. 4 The upper curve shows the periodogram of the time 
series {£"*(£;)}. The periodogram values are presented on 
the left axis. The lower curve illustrates the autocorrelation 
function of the same time series. The autocorrelation values 
are shown in the right axis. 



3m £ Ic, ds~^{m) = maxds^(Z) 



2M 



eC m . (8) 



To verify the implication (8) firstly it is necessary to 
evaluate the sets C; for I £ Ic of the values of r for 
which the autocorrelation function and the cosine func- 



tion with the period 



2M 



are positive and the sets 



Bj ,j £ Ib of values of r for which the autocorrelation 



2iU 



function and the cosine function with the period 
arc negative. Secondly, a percentage of the contribu- 
tion of the sum of products of positive values of c T and 
cos(^ : ) to the sum of positive products of the values 
of c T and cos(^ : ) should be evaluated. As a result the 
indexes ds^(l) for each set Ci where I is the index from 
the set Ic are obtained. Thirdly, from all sets Ci such 
that / £ Ic the set C m for which the index ds^(l) is 
the greatest should be chosen. 

The implication (8) is true when the set C m in- 

This means that 



2M 



eludes the considered period 
the greatest contribution of positive values of the au- 
tocorrelation function and positive cosines with the pe- 
riod 



2,1/ 



to the periodogram value Sk is caused by 

the sum of positive products of (X; — X)(Xi+ T — X) for 
each r6([f]-|, m + 1). 

When the implication (8) is false, the peak Sk is 
mainly created by the sum of positive products of {Xi — 
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Fig. 5 The upper curve shows the periodogram of the 
time series {S d (ti) — S d (ti)}. The periodogram values are 
presented on the left axis. The lower curve illustrates the 
autocorrelation function of the same time series (solid line). 
The dotted lines represent two standard errors of the auto- 
correlation function. The dashed horizontal line shows the 
zero level. The autocorrelation values are shown in the right 
axis. 



X)(X i+T - X) for each r 6 
where n is a multiple or a divisor of k. 

It is necessary to add, that the DE method should 
be applied to the periodograms peaks, which probably 
exist because of the echo-effect. It enables one to find 
such parts of the autocorrelation function, which have 
the significant contribution to the considered peak. The 
fact, that the conditions (7) and (8) are satisfied, can 
unambiguously decide about the existence of the con- 
sidered periodicity in the given time series, but if at 
least one of them is not satisfied, one can doubt about 
the existence of the considered periodicity. Thus, in 
such cases the sentence 'the peak can not be treated as 
true' should be used. 

Using the DE method it is necessary to remember 
about the power of the set /. If M is too large, errors 
of an autocorrelation function estimation appear. They 
are caused by the finite length of the given time series 
and result additional peaks of the periodogram 

occur. If M is too small, there are less peaks because 
of a low resolution of the periodogram. In applications 
M < ^ is used. In order to evaluate the value M 
the FFT method is used. The periodograms computed 
by the BT and the FFT method are compared. The 
conformity of them enables one to obtain the value M. 



4 Data analysis 

In this paper the sunspot activity data (August 1923 - 
October 1933) provided by the Greenwich Photohelio- 
graphic Results (GPR) arc analysed. Firstly, I consider 
the monthly sunspot number data. To eliminate the 11- 
year trend from these data, the consecutively smoothed 
monthly sunspot number (R(ti)) is subtracted from the 
monthly sunspot number (R(ti)) where the consecutive 
mean R(ti) is given by 



= i5 E for 



,N. 
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The values R(U) for i = 1,...,6 and i = N — 
6, . . . , N arc calculated using additional data from last 
six months of cycle 15 and first six months of cycle 17. 

Because of the north-south asymm etry of various so- 
lar indices ( Vizoso fc Ballester 1990l) . the sunspot ac- 
tivity is considered for each solar hemisphere separately. 
Analogously to the monthly sunspot numbers, the time 
scries of sunspot areas in the northern and southern 
hemispheres with the spacing interval At = 1 rotation 
are denoted. In order to find periodicities, the following 
time series are used: 
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{S n (ti) - S n (U)} 



{S S (U) - S S (U)} 



{R{U) - R(u)} 



- sunspot area fluctuations of the one 
rotation time interval in the north- 
ern hemisphere (N = 136), 

- sunspot area fluctuations of the one 
rotation time interval in the south- 
ern hemisphere (TV = 136), 

monthly sunspot number fluctua- 
tions (N = 122). 

In the lower part of Figure [1] the autocorrelation 
function of the time series for the northern hemisphere 
{S n (ti) — Sn(U)} is shown. It is easy to notice that 
the prominent peak falls at 17 rotations interval (459 
days) and c T for r € [3, 6] rotations ([81, 162] days) are 
significantly negative. The pcriodogram of the time se- 
ries {S n (U) — S n (ti)} (see the upper curve in Figures 
[T]) does not show the significant peaks at r = 5, 6 rota- 
tions (135, 162 days), but there is the significant peak 
at r = 9 (243 days). The peaks at t = 9, 17 are close 
to the peaks of the autocorrelation function. Thus, the 
result obtained for the periodicity at about 155 days are 
contradict to the results obtained for the time series of 



daily sunspot areas (jCarbonell k, Ballesterlll992l ). 

For the southern hemisphere (the lower curve in Fig- 
ure ^ c T for r € [2,7] rotations ([54, 189] days) is not 
positive except t = 5 (135 days) for which C5 = 0.02 
is not statistically significant. The upper curve in 
Figures [5] presents the periodogram of the time series 
{S S (U) — S S (U)}. This time series does not contain a 
Markov-type persistence. Moreover, the Kolmogorov- 
Smirnov test and the Fisher test do not reject a null 
hypothesis that the time series is a white noise only. 
This means that the time series do not contain an 
added deterministic periodic component of unspecified 
frequency. 

The autocorrelation function of the time series 
{R(U) — R(ti)} (the lower curve in Figure [3]) has only 
one statistically significant peak for r = 16 months 
(480 days) and negative values for r 6 [3,13] months 
([90, 390] days). However, the periodogram of this time 
series (the upper curve in Figure [3]) has two significant 
peaks the first at 15.2 and the second at 5.3 months 
(456, 159 days). Thus, the periodogram contains the 
significant peak, although the autocorrelation function 
has the negative value at r = 5 months. 

To explain these problems two following time series 
of daily sunspot areas are considered: 
{S d (ti)} - daily sunspot areas in the whole solar disk 
(N = 3653), 

{S (U) — S d (ti)} - daily sunspot area fluctuations, 
where 



i+182 



The values S d (U) for i = 1, . . . , 182 and i = N - 
182, ...,N are calculated using additional daily data 
from the solar cycles 15 and 17. 

The comparison of the functions c T of the time se- 
ries {S d (ti)} (the lower curve in Figured]) and {S d (ti} — 
S d (ti)} (the lower curve in Figure [5]) suggests that the 
positive values of the function c T of the time series 
{S d (ti)} in the interval of [61, 209] days could be caused 
by the 11-year cycle. This effect is not visible in the 
case of periodograms of the both time series computed 
using the FFT method (see the upper curves in Fig- 
ures S] and [5]) or the BT method (see the lower curve 
in Figure |6|) . Moreover, the periodogram of the time 
scries {S d (ti) — S d (U)} has the significant values at 
r = 117,158,203 days, but the autocorrelation func- 



tionis negative at these points. ICarbonell fc Ballester 
(|1992i) showed that the Lomb-Scargle periodograms fo r 
the both time series fsee lCarbonell &: Ballesterl (| 19921) . 
Figures 7 a-c) have a peak at 158.8 days which stands 
over the FAP level by a significant amount. Using the 
DE method the above discrepancies are obvious. To es- 
tablish the M value the periodograms computed by the 
FFT and the BT methods are shown in Figure |6] (the 
upper and the lower curve respectively). For M = 1000 
and for periods less than 166 days there is a good 
comformity of the both periodograms (but for periods 
greater than 166 days the points of the BT pcriodogram 
are not linked because the BT periodogram has much 
worse resolution than the FFT periodogram (no one 
know how to do it)). For M = 1000 and r = 154 the 



value of k is 13 ( 



2M 



153). 



satisfied because 2ds^ — 100 = 26% 



The inequality (7) is 

This means that 
the value of S13 is mainly created by positive values 
of the autocorrelation function. The implication (8) 
needs an evaluation of the greatest value of the index 

e solar data contain the most 
27 days because of the so- 



Thus, although ds~£ (0) = max; e / c ds~£(l) 



(9) 



ds~£(l) where I G Ic, but the solar data contain the most 
prominent period for r 
lar rotation. 

for each k < 13, all sets Ci (see (5) and (6)) with- 
out the set Co (see (4)), which contains r G [0,38], 
arc considered. This situation is presented in Figure 
[7] In this figure two curves y = cos^^- and y = c T 
are plotted. The vertical dotted lines evaluate the in- 
tervals where the sets Ai, Bj, Ci (for i, j, I = 1, . . . , k) 
are searched. For such C/ two numbers are written: 
in parentheses the value of dst(l) for the time series 
{S d (ti)} and above it the value of ds£(l) for the time 
series {S d (ti) — S d (U)}. To make this figure clear the 
curves are plotted for the set {r : r = 1, . . . , 500} only. 
(In the right bottom corner information about the val- 
ues of dst(l) for the time series {S d (ti) — S d (ti)}, for 
t = 501, . . . , 1000 arc written.) The implication (8) is 
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period (days) 

Fig. 6 The upper curve shows the FFT periodogram of the 
time series {S (U) — S d (ti)}. The FFT periodogram values 
are presented on the left axis. The lower curve illustrates 
the BT periodogram of the same time series (solid line and 
large black circles). The BT periodogram values are shown 
in the right axis. 




Fig. 7 Autocorrelation function of the time series {S (U) — 
S d (ti)} and the cosine function for k = 13 (the period at 
about 154 days). The horizontal line (dotted line) shows 
the zero level. The vertical dotted lines evaluate the in- 
tervals where the sets Ai, Bj, Ci (for i,j,l = 1, . . . , fe) are 
searched. The percentage values show the index ds^(l) for 
each Ci for the time series {S d {ti) — S d (ti)} (in parentheses 
for the time series {S d (ti)}). In the right bottom corner 
the values of ds~£(l) for the time series {S (tj) — S d (ti) }, for 
t = 501, . . . , 1000 are written. 



not true, because ds^{m) = 14.8% for m = 6. There- 
fore, [2M] = 153 i C 6 = [423,500]. Moreover, the 
autocorrelation function for r £ [115, 192] is negative 
and the set C2 is empty. Thus, ds^(2) = 0%. On 
the basis of these information one can state, that the 
periodogram peak at r w 154 days of the time se- 
ries {S d (U) — S d (ti)} exists because of positive c T , but 
for t from the intervals which do not contain this pe- 
riod. Looking at the values of ds^(l) of the time series 
{S d (ti)}, one can notice that they decrease when r in- 
creases until r = 346. This indicates, that when t 
increases, the contribution of the 11-year cycle to the 
peaks of the periodogram decreases. An increase of the 
value of ds^(l) is for I = 6 for the both time series, al- 
though the contribution of the 11-year cycle for the time 
series {S d (ti)} is insignificant. Thus, this part of the 

autocorrelation function (dst(6) = max; = i s dst(l) 

for the time series {S d (t t ) - S d (U)}) influences the fc-th 
peak of the periodogram. This suggests that the period- 
icity at about 155 days is a harmonic of the periodicity 
from the interval of [400, 500] days. 

The described reasoning can be carried out for other 
values of the periodogram. For example, the condition 
(8) is not satisfied for ft = 8, 9, 10 (250, 222, 200 days). 
Moreover, the autocorrelation function at these points 
is negative. These suggest that there are not a true peri- 
odicity in the interval of [200, 250] days. It is difficult to 
decide about the existence of the periodicities for ft = 6 
(333 days) and ft = 7 (286 days) on the basis of above 
analysis. The implication (8) is not satisfied for ft = 7 
and the condition (7) is not satisfied for ft = 6, although 
the function c T of the time series {S (U) — S d (ti)} is 
significantly positive for r = 286, 333. The conditions 
(7) and (8) are satisfied for k = 4 (Figure[S]) and k = 5. 
Therefore, it is possible to exist the periodicity from 
the interval of [400, 500] days. Similar results were 




Fig. 8 Same as in Figure 7, but for k = 4 (the 500-day 
period ) 
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925 950 975 1000 1025 1050 1075 

time (rotations) 

Fig. 9 Two upper curves show sunspot areas of the one 
rotation time interval in the whole solar disk S(ti) (solid 
line) and consecutively smoothed sunspot areas of the one 
rotation time interval S(U) (dotted line). Both indexes are 
presented on the left axis. The lower curve illustrates fluc- 
tuations of the sunspot areas S(U) — S(tC). The dotted and 
dashed horizontal lines represent levels zero and p respec- 
tively. The fluctuations are shown on the right axis. 



also obtained by iLeanl (|l990l ) for daily sunspot num- 
bers and daily sunspot areas. She considered the means 
of three periodograms of these indexes for data from 
N = 31 years and found statisticall y sign i ficant peaks 



from t he interval of [400 500] ( see |Lean| (jl99fj) . Fig 



ure 2). iKrivova fc Solankil (|2002l ) studied sunspot areas 
from 1876-1999 and sunspot numbers from 1749-2001 
with the help of the wavelet transform. They pointed 
out that the 154-158-day period could be the third har- 
monic of the 1.3-year (475-day) period. Moreover, the 
both periods fluctuate considerably with time, being 
stronger during stronger sunspot cycles. Therefore, the 
wavelet analysis suggests a common origin of the both 
periodicities. This conclusion confirms the DE method 
result which indicates that the periodogram peak at 
r = 154 days is an alias of the periodicity from the 
interval of [400,500]. 



5 The periodicity at about 155 days during the 
maximum activity period 

In order to verify the existence of the periodicity at 
about 155 days I consider the following time series: 
{S(ti) — S(ti)} - sunspot area fluctuations of the one 
rotation time interval in the whole 
solar disk (N = 136), 
{S?i(ti) — S%(ti)} - daily sunspot area fluctuations in 
the northern hemisphere from the 
maximum activity period (Jan- 
uary 1925 - December 1930, N = 

1478), 

{S*(U) ~ Sf(U)} - daily sunspot area fluctuations in 
the southern hemisphere from the 
maximum activity period (N = 
1478). 



The value S(U) is calculated analogously to R(U) 
(see Sect. 4). The values S^{ti) and Sf(U) are evalu- 
ated from the formula (9). In the upper part of Figure 
[9] the time series of sunspot areas (S(U)) of the one 
rotation time interval from the whole solar disk and 
the time series of consecutively smoothed sunspot ar- 
-6 S(U)) are showed. In the lower 



eas(5( ii ) = iES 



part of Figure [9] the time series of sunspot area fluc- 
tuations {S(ti) — S(ti)} is presented. On the basis of 
these data the maximum activity period of cycle 16 is 
evaluated. It is an interval between two strongest fluc- 
tuations e.a. J = [966, 1019] rotations. The length of 
the time interval J is 54 rotations. If the about 155- 
day (6 solar rotations) periodicity existed in this time 
interval and it was characteristic for strong fluctuations 
from this time interval, 10 local maxima in the set of 
Fi = S(ti) — S(ti) would be seen. Then it should be 
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necessary to find such a value of p for which Fi > p for 
i G J and the number of the local maxima of these val- 
ues is 10. As it can be seen in the lower part of Figure |9] 
this is for the case of p = 150 (in this figure the dashed 
horizontal line is the level of F t = 150). Figure [TU] 
presents nine time distances among the successive fluc- 
tuation local maxima and the horizontal line represents 
the 6-rotation periodicity. It is immediately apparent 
that the dispersion of these points is 10 and it is diffi- 
cult to find even few points which oscillate around the 
value of 6. Such an analysis was carried out for smaller 
and larger p and the results were similar. Therefore, 
the fact, that the about 155-day periodicity exists in 
the time series of sunspot area fluctuations during the 
maximum activity period is questionable. 

To verify again the existence of the about 155-day 
periodicity during the maximum activity period in each 
solar hemisphere separately, the time series {S n (ti) — 
S n (U)} and {S S {U) — S S (U)} were also cut down to 
the maximum activity period (January 1925-December 
1930). The comparison of the autocorrelation functions 
of these time series with the appriopriate autocorrela- 
tion functions of the time series {S n (U) — S n (U)} and 
{S s (ti) —S s (ti)}, which are computed for the whole 11- 
year cycle (the lower curves of Figures [1] and [5]) , indi- 
cates that there are not significant differences between 
them especially for r=5 and 6 rotations (135 and 162 
days)). This conclusion is confirmed by the analysis of 
the time series {S^iU) — S%(U)} for the maximum ac- 
tivity period. The autocorrelation function (the lower 
curve of Figure [TTj) is negative for the interval of [57, 
173] days, but the resolution of the periodogram is too 
low to find the significant peak at r = 158 days. The 
autocorrelation function gives the same result as for 
daily sunspot area fluctuations from the whole solar 




- -0.2 



200 250 300 
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Fig. 11 Same as in Fig. 5, but for the time series {S%(ti) - 
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Fig. 10 Nine time distances among the successive local 
maxima of the values of Fi > p. The horizontal line repre- 
sents the 6-rotation (162-day) period. 



Fig. 12 Same as in Fig. 5, but for the time series {Sf (ti) - 
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disk ({S d (ti) — S d (ti)}) (see also the lower curve of Fig- 
ures [5]). In the case of the time series {S s (ti) — S s (ti)} 
C5 is zero for the fluctuations from the whole solar cy- 
cle and it is almost zero (c T = 0.02) for the fluctuations 
from the maximum activity period. The value cq is neg- 
ative. Similarly to the case of the northern hemisphere 
the autocorrelation function and the periodogram of 
southern hemisphere daily sunspot area fluctuations 
from the maximum activity period {S~ (ij) — Sf{ti)} are 
computed (see Figure fT2|) . The autocorrelation func- 
tion has the statistically significant positive peak in the 
interval of [155, 165] days, but the periodogram has too 
low resolution to decide about the possible periodicities. 
The correlative analysis indicates that there arc posi- 
tive fluctuations with time distances about 155 days in 
the maximum activity period. 

The results of the analyses of the time series of 
sunspot area fluctuations from the maximum activ 



ity pe riod are contradict with the conclusions of iLean 
(1990). She uses the power spectrum analysis only. 



The periodogram of daily sunspot fluctuations contains 
peaks, which could be harmonics or subharmonics of 
the true periodicities. They could be treated as real 
periodicities. This effect is not visible for sunspot data 
of the one rotation time interval, but averaging could 
lose true periodicities. This is observed for data from 
the southern hemisphere. There is the about 155-day 
peak in the autocorrelation function of daily fluctua- 
tions, but the correlation for data of the one rotation 
interval is almost zero or negative at the points r = 5 
and 6 rotations. Thus, it is reasonable to research both 
time series together using the correlative and the power 
spectrum analyses. 



6 Conclusion 



The following results are obtained: 

(1) . A new method of the detection of statistically sig- 
nificant peaks of the periodograms enables one to 
identify aliases in the periodogram. 

(2) . Two effects cause the existence of the peak of the 
periodogram of the time series of sunspot area fluc- 
tuations at about 155 days: the first is caused by 
the 27-day periodicity, which probably creates the 
162-day periodicity (it is a subharmonic frequency 
of the 27-day periodicity) and the second is caused 
by statistically significant positive values of the au- 
tocorrelation function from the intervals of [400, 500] 
and [501, 1000] days. 



(3) . The existence of the periodicity of about 155 days 
of the time series of sunspot area fluctuations and 
sunspot area fluctuations from the northern hemi- 
sphere during the maximum activity period is ques- 
tionable. 

(4) . The autocorrelation analysis of the time series of 
sunspot area fluctuations from the southern hemi- 
sphere indicates that the periodicity of about 155 
days exists during the maximum activity period. 
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